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FOURIER SERIES OPERATING PACKAGE 


by 

Milton L. Charnow 
Goddard Space Flight Center 


INTRODUCTION 


The Hansen Satellite Theory as modified by Musen 1 involves various manipulations 
of "literal" Fourier series, as such, before final numerical evaluation. To program this 
formulation it is necessary to represent cosines, sines, and constants in a special man- 
ner and thus manipulate the various Fourier series. 


SPECIAL REPRESENTATION OF A SINE, COSINE, 
OR CONSTANT 


It was decided to represent any term, including the constant term, of the Fourier 
series we are concerned with by two consecutive eight-place floating-point words. The 
terms of this series have the general form 

A n cos (iF + jE + kW + 1U) 
or 

A n sin (iF + jE + kW + 1U) 

and the values of F, E, W, and U are not used until final numerical evaluation. The first 
eight-place floating-point word represents the coefficient A n in an entirely standard man- 
ner. The second eight-place floating-point word represents the sine or cosine term of up 
to four arguments in an artificial manner. 


1 Mus en, P., "A Modified Hansen's Theory as Applied to the Motion of Artificial Satel- 
lites, " NASA Technical Note D-492, November I960. 
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Consider a term such as A n cos (iF + jE + kW + 1U). U itil this term is evaluated we 
are concerned only with the values of A n , i, j, k, and 1. The coefficient A n is represented 
by the first of the two eight-digit floating-point words. The Lour coefficients of the argu- 
ments (i, j, k, and 1) are each represented by two digits of the second eight-place floating- 
point word. The first two digits represent i and are normalized to 00; the next three pairs 
of digits represent j, k, and 1, and are normalized to 50. Thus, the value of i may vary 
from 0 to 99, while j, k, and 1 may vary from -49 to +49. A cosine is denoted by a plus 
sign; a sine by a minus sign, that is, 1 cos (OF + 0E + 0W + 0U) or cos (0) becomes 
+ 10000000 + 01, +00505050 + 08. Any constant term can be represented as A n cos (0). A 
few additional examples will be helpful (see also Appendix A): 


Conventional 


Spec ial 


1/4 cos (IF + 2E - 3W - 2U) 
1/4 sin (OF + IE - 0W + 2U) 
1/4 
0 


25000000 + 00, 
25000000 + 00, 
25000000 + 00, 
00000000 + 00, 


+01524748 + 08 
-00515052 + 08 
+00505050 + 08 
+00505050 + 08 


The convention of having the first non-zero coefficient of the argument word positive was 
adopted. Cos(-x) becomes cos(x) and sin(-x) becomes -sin(x). Examples are: 

(1/2) cos (OF - 3E + 1W + 0U) becomes (1 2) cos (OF + 3E - 1W + 0U) 
and the special representation is 50000000 + 00, + 00534950 + 08 

(1/2) sin (OF - 3E + 1W + 0U) becomes -(1/2) sin (OF + 3E - 1W + 0U) 
and the special representation is -50000000 + 00 - 00534950 + 08. 

The first location address of a series contains the num ier of terms of the series. 

A series of n terms would be represented by 2n + 1 words the first of which would be 
the number n. 


THE FOURIER OPERATING PACKAGE 

The following series operations are performed by the Fourier Operating Package: 

• Multiplication 

• Addition and Subtraction 

• Differentiation 

• Integration 

• Bar (Special operation used in the Hansen satellite theory) 

• Scalar Multiplication 

• Coefficient Extraction 

• Series Evaluation 


Q h 
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Multiplication 

Multiplication of two series, where the terms are of the general form described 
earlier and the values of F, E, W, and U are not used until the final numerical evalua- 
tion, is according to the conventional trigonometric identities: 

A cos X-B cos Y = (AB/2) cos (X+Y) + (AB/2) cos (X-Y) 

A cos X-B sin Y = (AB/2) sin (X+Y) - (AB/2) sin (X-Y) 

A sin X-B cos Y = (AB/2) sin (X+Y) + (AB/2) sin (X-Y) 

A sin X-B sin Y = (AB/2) cos (X+Y) + (AB/2) cos (X-Y). 

The Fourier Multiplication routine is composed of three major sections: the multiplier, 
the collapser, and the arranger. 

The Multiplier 

Two Fourier series such as 

AjAj + A 2 A^ + AgAg + . . . + v£ (SerieS A) 

aild ^ ^ ^ 

B 1 B 1 + B 2 B 2 + B 3 B 3 + • • ' + B m B m (Series B) 

which are to be multiplied are arranged in descending order of the absolute values of the 
coefficients, that is, 

Af > A 2 > > . . . > A n and Bj > B 2 > B 3 > . . . > B m 

To facilitate further discussion, we shall denote any term in the A series A X A£, 

* * 
any term in the B series as B y B y , and any term inthe resultant series by C z C z , where 

A , B and C are the coefficient words and A , B , and C z are the argument words. 
a y z y 

The multiplication of the A series by the B series proceeds as follows: The first 
term in the A series is multiplied by each term in the B series, then the second term in 
the A series is multiplied by each term in the B series, and so on until each term in the 
A series has been multiplied by each term in the B series. For example^] A x By| is compared 
with some numerical criterion e. If | A x By| > e, then C Z C Z and C( z + 1 ) C( z + 1 ) are gen- 
erated according to the trigonometric formulas already stated. If | A x By| < e, then the 
value of y is examined. If y > 1 (i.e., B y is any term other than the first term), A x is 
replaced by A^ x + ^ and A^ x + is multiplied by Bj. If y = 1, the multiplication of the 
two series is terminated since any further A^ x + ^ By will also be less than e. The mul- 
tiplication process continues until A^A^ has been multiplied by B m B ^ 1 unless the 
numerical criterion or space limitations intervene. 
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The CoUapser 

Every multiplication generates two terms of two words each.* The purpose of the 
collapser is to combine like argument terms. Each argument term C z is compared with 
each other argument term previously generated and stored. If C z equals any other 
argument term, the corresponding coefficient terms are added. Thus, there is no dupli- 
cation of terms. 

The Arranger^ 

The final step in the multiplication is the arranging of the terms of the series. 

| Ci | is compared with | C2 1, | C3 1 , etc. If | C z | > | Cj |, then Cj is replaced by C z and 
C| is replaced by C*. The process continues until the terms are arranged, in descending 
order, according to the absolute value of the coefficients. 


Addition and Subtraction 


Addition or subtraction of two Fourier series is primarily a process of com- 
paring argument terms and adding the coefficients of like terms. Aj^ is compared suc- 
cessively with b| through A*) with B| through B ri , etc., until A n has been compared 
with B^. If A* = B*. the sum of the coefficients (A x 1- B y ) and the argument term A* are 
stored, and B y and B y are replaced by zeros If A^ does not equal any B y , both A x and 
A* are stored. After all terms in the A series have been compared with all terms in the 
B series, the remaining B series terms are stored. 

Subtraction is accomplished in like manner after changing the signs of each coeffi- 
cient term in the B series. The resultant series in each case is processed through the 
arranger. 


Differentiation 


Differentiation, in this application, is with respect to the F variable. Thus, 

(lllr) a sin < iF + i E + kW + 1U ) = iA cos < iF + i E + kW + iu). 


Example : 


(J^j sin (3F + 2E - 3W + U) = +3 ccs (3F + 2E - 3W + U) 
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and 

[10000000 + 01, - 03524751 + 08 becomes +30000000 + 01, +03524751 + 08. 

On completion of the differentiation, the resultant series is processed thru the 
arranger. 


Q 


Integration 

Integration, in this application, is with respect to E. However, W is also a func- 
tion of E. Thus, 

Ja cos (iF + jE + kW + lU)dE = — j + A C2 ^ sin (iF + jE + kW + 1U). 

Example : 

With ci = 1 and c 2 = 1, 

J[30000000 + 01, + 01525250 + 08]dE becomes 75000000 + 00, -01525250 + 08. 

The resultant integrated series is also processed through the arranger. 


Bar Operation 

The bar operation is a special function in the Hansen Satellite Theory. It con- 
sists of adding the coefficient of the F argument to the coefficient of the E argument and 
substituting zero for the F coefficient. Thus, 

A cos (iF + jE + kW + 1U) after barring becomes A cos (OF + (i + j) E + kW + 1U). 

Example : 

50000000 + 00, + 02534850 + 08 after barring becomes 50000000 + 00, + 00554850 + 08. 


Scalar Multiplication 

Scalar multiplication is the multiplication of the coefficient A n of each term by 
a constant. 
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Coefficient Extraction 

In the Hansen Satellite Theory it is sometimes necessary to use the coefficient 
of some term of a Fourier series such as a sine IF term, cosine 2E term, or the con- 
stant term of a series. Let us assume it is necessary to use the constant term. If there 
is a constant term in that series, it will be the multiplier of the cos (0), or in special 
form, that A n word which multiplies 00505050 + 08. We successively compare each 
argument term in the series with cos (0) and extract that A n which multiplies the argu- 
ment term cos (0). If no argument term of the series is cos (0), a normalized zero 
(10000000 + 00 + 00000000 + 00 + 00505050 + 08), is printed. 


Series Evaluation 

The numerical values F, E, W, and U are only employed in the Series Evalua- 
tion Routine. 

To evaluate a Fourier series, the numerical values of i, j, k, and 1 are multi- 
plied by the numerical values of F, E, W, and U, and the sum iF + jE + kW + 1U is de- 
termined. The sine or cosine of iF + jE + hW + 1U is multiplied by the coefficient A and 
the terms are added. 


CONCLUDING REMARKS 

Appendix A is the special representation of a nine term series. Appendix B presents 
flow charts and Appendix C a listing of instructions for the program packages for multi- 
plication (including the collapser and arranger), addition or subtraction, differentiation, 
integration, bar operation, scalar multiplication, coefficient extraction, and the evalua- 
tion of the final series. Because this program was writter in Mystic Code for the IBM 
709, an explanation of Mystic Code is given in Appendix D. 

The Fourier operating package can be used with any theory that involves representa- 
tions of functions by Fourier series. It can also be modified to operate with polynomials 
of the form 


8101 - 



7 


ACKNOWLEDGMENTS 

The author is indebted to Messrs. R. G. Kelly and T. P. Gorman for their aid in the 
construction of the package and for the translation into Mystic Code and to Aileen Mario 
for preparing the flow charts. 




-1078 


9 


Appendix A 


Nine Term Series in Standard and Special Form 


Q The following is a nine term series presented in standard and in special form. Note 

that in the special form, the first location (address) contains the number of terms in the 
series. 


Series in standard representation 

.29467121 cos (0) 

+.00010496334 cos (OF + 0E + 2W + 0U) 
+.00005252596 sin (OF + IE - 1W + 0U) 
+.000019845618 cos (OF + IE - 2W + 0U) 
+.0000066329604 cos (OF + IE + 2W + 0U) 

+ .0000020107054 cos (OF + 2E - 2W + 0U) 
-.00000036004597 sin (OF + 0E + 1W + 0U) 
-.000000055052357 sin (OF + OE + 3W + OU) 
-.000000031090653 sin (OF + IE + 1W + OU) 


Series in special representation 


90000000 + 01 

+29467127+00, +00505050+08 
+ 10496334-03, +00505250+08 
+52525962-04, -00514950+08 
+19845618-04, +00514850+08 
+66329604-05, +00515250+08 
-20107054-05, +00524850+06 
-36004597-06, -00505150+08 
-55052357-07, -00505350+08 
-31090653-07, -00515150+08 




Appendix B 


Flow Charts 

The following are the flow charts for the op- 
erating packages for multiplication (including 
the collapser and the arranger), addition or 
subtraction, bar operation, differentiation, in- 
tegration, and series evaluation. 
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Flow Chart for Multiplication (Continued) 
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Flow Chart for Multiplication (Continued) 
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Appendix C 


Listing of Instructions 


The following is a listing of instructions for 
multiplication (including the collapser and the 
arranger), addition or subtraction, differentia- 
tion, integration, bar operation, scalar multi- 
plication, and coefficient extraction. 

The square root instructions are for stand- 
ard number representations and therefore have 
not been discussed in the report. 
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Appendix D 

Mystic Code 


The Fourier Operating Package described in 
this report was written for the IBM 709 in Mys- 
tic Code. The following is an explanation of the 
Mystic Code. 
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